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ABSTRACT 


An  algorithm  for  predicting  the  Max i mum  Usable  Frequency 
between  two  points  on  the  surface  of  the  earth  for  High 
Freauency  communications  is  implemented  on  a  TI-59  hand-held 
programmable  calculator* 


1*  INTRODUCTION 

A  n  a  I  a  o  r  i  t  h  in  for  ppedictin  s*  t  h  e  M  a  x  i  m  u  m  U  s  3  b  1  e 
Freauency  (MUF)  in  the  High  Frequency  radio  band  has  been 
developed  at  the  Naval  Ocean  Systems  Center*  Sari  Diego? 
California*  by  R*B*  Rose  3nd  J*N*  Martin  C 1 1 ?  and  has  been 
implemented  in  both  FORTRAN  IV  and  BASIC  computer  langusges 
under  the  name  MINIMUF-3»5«  The  algorithm  has  been 
incorporated  into  the  CLASSIC  PROPHET  propagation  analysis 
system*  and  has  been  demonstrated  to  provide  a 
f ield-deployable  capability  for  computing  HF  propagation 
paths  with  micro-computer  resources* 

2.  OBJECTIVE 

The  essential  portions  of  MINIMUF-3*5  are  contained  in 
approximately  80  BASIC  language  statements  *  reproduced  in 
appendix  A*  which  suggests  a  potential  for  being  encoded  in 
hand-held  programmable  calculators*  H*F*  Hite  of  Hughes 
Aircraft  Company  adapted  a  restricted  version  of  MINIMUF--3  ♦  0 
for  the  HP-67  calculator  that  computes  MUFs  for  1-hop 
F- layer  propagation  and  is  thus  limited  to  station 
separations  of  4000  kilometers  or  less*  A  complete 
implementation  of  MINIMUF*  capable  of  2-hop  path  prediction* 
W3s  desired* 


3*  PROGRAM  DEVELOPMENT 

A  TI-59  calculator  program  was  adapted  directly  from 
the  Rose  and  Martin  BASIC  program  listing  with  some  logic 
changes  required  to  accomodate  the  memory  constraints  of  the 
calculator*  The  program  listing*  summary  of  logic  changes* 
operating  instructions r  and  test  C3se  are  contained  in 
appendices  B>  C?  Br  and  E?  respectively*  The  program  listing 
is  annotated  with  cross-references  to  statements  in  the 
BASIC  listing  for  assistance  in  tracing  logic  flow*  The  hulk 
of  the  development  task  was  to  minimize  program  and  data 
storage  reaui rements ♦  Several  iterations  finally  produced 
the  current  version  which  reaui res  800  program  steps  and  20 
data  registers f  the  full  capacity  of  the  normally-accessible 
calculator  attributes* 

4*  PERFORMANCE 

The  TI-59  program  was  optimized  for  storage  r  that  being 
the  overriding  constraint*  Execution  time  for  single-hop 
predictions  is  approximately  50  seconds*  Were  the  program 
optimized  for  speed?  given  a  larger  storage  capacity? 
execution  time  could  conceivably  approach  40  seconds*  The 
test  case  contained  in  appendix  E  is  the  identical  test  case 
promulgated  with  the  Rose  and  Martin  report y  and  yields  the 
same  results  on  the  TI-59*  Notice ?  however?  that  the  MUF 
limit  of  32  MHz  has  been  raised  to  50  MHz*  thereby  providing 


for  HF  path  prediction  during  periods  of  high  solar  flux 
when  customary  band  limitations  are  exceeded.  The  test  case 
should  be  executed  after  initially  programming  the  TI-59  by 
keystroke  to  ensure  correct  program  entry,  A  copy  of  the 
program  is  available  on  magnetic  cards  by  sending  two  blank 
cards  to  the  authors r  however r  the  user  is  cautioned  that 
magnetic  cards  are  not  guaranteed  to  be  transportable  among 
all  TI-59  calculators. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  capability  for  predicting  Maximum  Usable  Frequencies 
in  the  HF  radio  b3nd  can  be  realized  with  hand-held 
calculators r  and  therefore  can  be  widely  deployed  in  the 
field*  An  even  larger  capability  can  be  realized  by 
utilizing  custom  made  modules  for  the  TI-59.  Each  module 
provides  for  5000  program  steps  in  addition  to  the  normal 
960  maximum  step  capacity  of  the  calculator »  and  could 
provide  for  additional  portions  of  the  CLASSIC  PROPHET 
System  to  be  included  in  calculator  implementations r  such  as 
the  D-reSion  absorption  model.  Appendix  F  lists  several 
persons  to  contact  in  regard  to  having  custom  modules 
manufactured  for  the  TI-59. 


APPENDIX  A 


BASIC  MINIMUF-3.5  PROGRAM  LISTING 


The  BASIC  language  listing  of  MINIMUF-3.5  is  reprinted  here 
with  the  permission  of  R.B.  Rose*  and  serves  as  a  Suide  for 
the  TI-59  keystroke  listing. 


MIN1MUF-3.5  PKOC.K AM 


ieee  rem 

1818  K7-SIN<L1>*SIN<L2>+C0S<L1)*C0S(L2>*CQS<W2-M1> 
1828  G1«ACS<K7  MAX  -1+1.8E-3  MIN  1-1.8E-8) 
1938  K6-1.391G1 

1848  K6-K6  MAX  1 

1858  K5«l/K6 

1868  J9-188 

1878  FOR  Kl-l/<2tK6)  TO  1-1/<2*K6>  STEP  8.9999-1/K6 

1888  IF  K3-1  THEN  1188 

1090  K5-8.5 

1188  P»SIN<L2> 

1119  Q«C0S<L2) 

1129  A«<SIH<L1)-P*C0S<G1>)/<Q*SIN<G1>> 

1138  3-G13X1 

1148  OP3C0S<B>*QtSIN<B)*A 

1158  D-<C0S<B>-C*P)/<Q*SQR<1-C*C>> 

1168  D«ACS<D  MAX  -1+1.8E-5  MIN  1-1.8E-8) 

1178  M0«M2+SGN<SIN<W1-W2)>*D 

1188  IF  W0«>8  THEN  1288 

1190  W9-W8+P1 

1288  IF  M8<P1  THEN  1228 

1218  M8»W8-P1 

1229  L9«P9-ACS<C  MAX  -1+1.0E-5  MIN  1-1.8E-8) 

1238  Y1*8.8172*<18*<M8-1>«38.4+D6> 

1240  Y2*8.489*CQS<Y1> 

1250  K8-3.823M8+12+8. 13*<SIN<Y1 >*1.2*SIN<2*Y1>> 

1260  K8»K8-12*<l+SGN<K8-24>XSGN<ABS<K8-24>> 

1278  IF  C0S<L8+Y2>>-8.26  THEN  1350 

1280  K9»8 

1298  G8«8 

1380  N9«2.5*G1*K5 

1310  M9-H9  MIN  P0 

1329  M9«SIN<M9> 

1339  M9«1+2.5*H9*SQR<M9> 

1340  GO  TO  1590 

1358   X9«<-0.26*SIN<Y2>*SINa0>>/'<COS<Y2>*COS<L0>  +  1.0E-3> 

1360  K9«12-ATN<K9/SQR<ABS<1-K9*K9>>>*7. 639437 

1 370  T-K8-K9/2+ 1 2% < 1 -SGN < K8-K9/2 > ) *SGN < ABS  <  K8-K9/2 ) ) 

1388  T4«K8+K9/2-12*a*SGN<X8+K9/2-24mSGN<ABS<K8+K9/2-24>> 


1390  C0=ABS<COS<L0+Y2>> 

1480  T9«9.7*C8t9.6 

1418  IF  T9>0.1  THEN  1438 

1420  T9=0.1 

1430  M9=2.5*GUK5 

1440  M9=M9  HIM  P0 

1450  H9-SINCH9) 

1460  M9«1+2.5*H9*SQR<M9> 

1478  IF  T4<T  THEN  1300 

1480  IF  <T5-T>*<T4-T5>>0  THEN  1310 

1490  GO  TO  1648 

1588  IF  <T5-T4>*<T-T3>>0  THEN  1640 

1518  T6=T5+12*(1+SGN<T-T5>>*SGN<ABS<T-T5>> 

1529  G9»P«<T6-T>/K9 

1538  G8»PI$T9/K9 

1548  U»<T-T6>/T9 

1558  G8=C8*<SIN<G9>+G8*<EXP<U)-C0S<G9>)>/<H-G8*G8> 

1560  G7«C9*<G8*<EXP<-K9/T9»1 ) >tEXP< <K9-24>/2>/< 1*G8*G8> 

1570   IF  G0->G7  THEN   1390 

1580  G0*G7 

1590  G2»<H>S9/250>*H9*SQR<6*58*SQR<G8>> 

1600  G2«G2*<l-0.UEXP<<K9-24>/3>> 

1610  G2»G2SU«K1-SGN<L1>*SGN<L2>>*0.1> 

1629  G2«G2$<1-0.U<1+SGN<A8S<SIN<L0>)-COS<L0>>>> 

1630  GO   TO    1700 

1640  T6*T5+12*< 1+SGN<T4-T5> >*SGN< ABS<T4-T5> ) 

1650  G8«PI*T9/K9 

1660  U*<T4-T6)/2 

1670  U1—K9/T9 

1689  G0«C9*<G8*<EXP<U1 >+l > >*EXP<U>/<1+G8*G8> 

1698  GO  TO  1590 

1700  IF  G2>J9  THEN  1720 

1710  J9-G2 

1720  NEXT  Kl 

1730  J9-J9  MAX  2  H IN  32 

1740  RETURN 
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APPENDIX  B 


TI-59  PROGRAM  LISTING 


The  keystroke  listing  of  the  TI-59  program  implementation  of 
MINIMUF-3»5  follows*  Segments  of  TI-59  code  are 
cross-referenced  to  the  BASIC  program  listing* 


~i  r 


60    DEG    1515720 


050 
05  1 


■■IT      1080,90 


1300,10,20,30 


!~i  C7  O 


U 


1030,^0-50 


1060 


i  i :-! ;: 


:■-.:■- 

s_i  y 

-i 

04  1 

1. 1  cr 

= 

0  4  2 

.  S  £_ 

■::.  1  T 

0  -L  3 

0  1 

1 

"■'  i 

S 13  R 

04  F: 

07 

u  r 

H4£ 


•--■  j     I  •••  X 

42    3TD 


u4'-;i     r 


i       .-■  i 


Ob y         1       il 


!  j  9  5  95  - 

0 96  42  S T D 

0  9  7  0  9  i  j  9 

09  8  01  1 

I"I9.9  '-Fi 


1070 


1070  -  cont'd 


10 


100       43    RCL 


J.    !_'  u_ 

I  0  3 

I  0  4 

105 

1  u£ 

107 

1  0  :: 
i  o  9 

1  10 

i  1 


U  f"  U  I 


7  cr 


1  5Q 


1140 


1       "      '~l 

I'd. 
1      -,    ~ 


!  .-;  i 


E-"  Er  j. 


1130 


RCL       1100,10,20 


1220 


1150 


1  :-:4 


I  !"•■' 
h  .  1 


•  :■  9    L-  lJ  i 
95       = 

I.-  tr 


1  9 

1  o . 


19' 

19! 

1    Q  • 


J 


11 


200       22    INV       1160,70,80,90,1200,10        250 


i  ■-' 


1 4  1       n  i 


i ! :.-' 


t-,. 


1230 


1240 


1250 


2^6       I..!] 
247       lin 


V  '-J  S 


297       65 


249       43    RCL 


299       4:3    RCL 


12 


"I       II I 


'01 
!  0  2 
:  n :"! 


Ijv 
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1590 


;  I  J  h: 

;  n  q 


T2M 


IV 


I  i. 


i  ;  i 


ill  I 


1270 


•',  f-ll 


1600 


39    CDS 


•!■■  '—  T 

s-2 

C  t  n 

'':  lL  r~- 

1    !    ".' 

:_!  jl 

.']    -T    *~l 

"?  "7 

be 

347 

04 

04 

b  0 
0  0 

-..  0 

3  4  3 

0     1280,90 

':"'  6 
QQ 


I  i  ! 


t'l  i 


1610 


13 


4  0  0 

■ifl  i 


01 


4.h,n       it 


'-r-  i_!  ■-* 


.-.    i"l  Cl 


TS20" 


1350,60 


439 


K 1 


1700 


4. 4 


T  -J 


448 

4  49 


1   1 


;z    1070,1720 


u  u 


!  j'-l 

!J  0 
9lr! 


1730,^-0 


■-t  ~ 


■J    v 


496 

■  ■  c. 

49? 

1 1 1 

4  93 

02 

4  9  9 

95 

14 


5 0 0      42    S T  D 


101 


J  u  ^ 


43    RCL     1370 


*=:«=;  n 


i  r 


err  ■: 

C  !=j 

< 

-      — 

cr  C  ■-: 

•-  !     -1 

'..'     !■    - 

•—  ■-   £_ 

-'    *— 

i    *. 

cr  cr  -i 

.     i  .     i  .  ": 

fj          l 

1 

_•!_ 


=m7 


D  1  U       £  z 

511       32 
5  1  2       U  1 


1390,1^00,10,20 


"1 1  a       b . 


I  ■"-! 


1380 


1530,1650 


~  .-1    1 


~    7i: 


3    RCL    1^70 


543 
*-!  ii  si 

5  4  6 


01      i 


hi  "> 


6?  a? 


c,  Q  ^ 

cr,  9  "•" 
5  9  9 


05 


b9  a:' 
i  n      1  n 


15 


b  U  Li 

s-    :  1 

i-  :";  :'" 


650 

652; 


b  i  O 

h  i  4 


■S  -t'CT 


RCL     1510 


1550 


,-.  ,  •;  ■-) 

0  bi 

i"i  H      i 

640 

t~!     1 

i 

64  1 

54 

'\ 

642 

6  5 

6  4  3 

■J  wl 

'A ;  T 

644 

T    ..•   T 

645 

■7  cr 
.j 

- 

646 

43 

RCL 

647 

04 

04 

64  S 

9  ~i 

— 

649 

65 

1520 


A  -i  ;_: 

1  14 

69  1 

O  !-" 

6  9  2 

0 ... 

6 i_i  3 

''  T  O 

h  ;=l  4 

00 

6  9  5 

.-■  cr 

6  9  0 

42! 

697 

02 

6  9  8 

^5 

6  9  9 

C]  O 

rC    RCL     T5S0" 


16 


70  1 

702 


)  •-.       r:  f  I 


i  U  i 

7  0  8 

7  0  9 
7'  1  0 

"7  1    1 

;'    ii 
i'    i  O 


Ill  1 


1680 


1570,80 


1640,60 


MIN/MAX 

SUBROUTINE 


47 
4  ;•; 

69    DP 

4':-i 

S  5       f 

ro     >:v   ;;:t 
99       92    PIN 
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APPENDIX  C 


PROGRAM  LOGIC  MODIFICATIONS 


The  following  progrsm  logic  changes  were  made  to  the  BASIC 
version  of  MINIMUF-3.5  to  minimize  program  storage 
reoui  rements I 

(1)  The  computation  of  M9  at  lines  1300-.1.330  and  1430-1460 
is  independent  of  the  main  loop  and  was  moved  to  the 
beginning  of  the  program  t  immediately  following  the 
calculation  of  Gl  3nd  K5* 

(2)  The  loop  index  computation  at  line  1070  was  separated 
into  3  computation  of  a  new  variable r  KHOP*  which  is  either 
1  or  2?  and  3  follow-on  computation  of  Kl  from  KHOP  as  the 
first  item  within  the  loop* 

(3)  The  computation  of  G3  3t  lines  1530  3nd  1650  is  common 
to  both  le^s  of  the  preceeding  test  for  ( T5-T4 M T-T5) *  snd 
has  been  moved  ahead  of  that  test. 

(4)  The  test  at  line  1470  serves  to  reverse  the  sense  of  the 
following  tests  at  lines  1480  and  1500.  This  logic  has  been 
combined  into  a  single  test  on  ( T5-T4 ) ( T-T5 )  with  3 
weighting  fsctor  of  1  or  -1  to  reverse  the  sense  of  that 
test. 


IS 


<5)  The  computation  at  line  1360  was  simplified  to  the 
srcsine  by  way  of  trigonometric  identity. 

(6)  The  MUF  limit  established  at  line  1730  was  changed  to  50 
vice  32  to  accomodate  high  solar  flux  densities* 

(7)  Throughout  the  program*  explicit  calculation  and  storage 
of  variables  that  are  used  only  once  in  following  statements 
has  been  eliminated  to  conserve  on  data  register 
reeui rements ♦ 
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APPENDIX  Ii 


PROGRAM  OPERATING  INSTRUCTIONS 


The  following  instructions  must  be  followed  to  operate  the 
TI-59  version  of  MINIMUF-3»5: 

(1)  Repartition  the  calculator  for  800  program  steps  and  20 
data  registers  by  pressing  2/0P/17* 

(2)  Lo3d  the  4  memory  banks  from  program  C3rds  <2>  or  by 
entering  program  keystrokes  with  the  calculator,  in  the  LRN 
mode  ♦ 

(3)  Enter  input  d3ta  as  follows! 

(a)  Transmitter  North  Latitude  in   decimal   degrees   in 
R12  (range  -90  degrees  to  90  degrees) 

(b)  Transmitter   West  Longti tude  in  decimal  degrees  in 
R13  (range  0  to  360  degrees) 

(c)  Receiver  North  latitude  in  decimal  degrees   in   Rl-4 
(range  -90  to  90  degrees) 

(d)  Receiver  West  Longti tude  in  decimal  degrees  in  R15 
(range  0  to  360  degrees) 


20 


(e)  Month  in  R16  ( rar^e  1  to  12) 

(f)  Day  in  R17  ( ranse  1  to  31) 

(g)  Sunspot  Number  in  R19  ( ran^e  -  positive  number) 

(4)  Enter  Time  in  X-register  (range  0  to  24  decimal  hours) 

(5)  Press  A 

(6)  If  a  printer  is  attached!'  the  input  time  is  echoed  r 
followed  by  the  answer  (MUF)  after  50  seconds  for  single-hop 
and  100  seconds  for  two-hop  predictions* 

(7)  The  answer  is  displayed  in  the  X- register  when  the 
program  halts* 

(8)  All  input  Quantities  remain  undisturbed  in  R12  thru  R.19 
(time  is  stored  by  the  program  in  R13).  Another  MUF  for  a 
different  time  can  be  computed  directly  by  repeating  steps 
(4)  and  (5>« 
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APPENDIX  E 


TEST  CASE 


The  test  case  that  follows  was  provided  in  the  NOSC  report 
on  MINIMUF-3»5»  The  actual  printer  listing  of  input  and 
program  output  is  included* 


n,  n  (Time) 
CONTENTS      REG  .-^  .-,  ^-p) 


?.  1  = 
1  F>.f,. 


U. 

'  u. 


1  2 

(LI) 

(wi) 

14 

(L2) 

15 
1  f-. 

(W2 

(MO) 

(D6) 

i.  •_■ 

1  c 

j.  -•• 

(S9) 

INPUT 
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:5.  0  (etc) 


.-  3    ■_• 


<  -. 


20.  i.l 


OUTPUT 


APPENDIX  F 
TI-59  CUSTOM  MODULES 


The  following  points  of  contact  3re  provided  for  the  design 
3nd  manufacture  of  TI-59  custom  modules*  Both  Texas 
Instruments  and  Horizons  Technology  (contracting  with  TI) 
provide  software r  emulator r  and  consulting  support  for  the 
design  and  production  of  custom  modules  ♦  Costs  Quoted  at 
this  time  are  approximately  $12*000  for  a  minimum  order  of 
250  modules.  This  includes  emulator  support  and  manufacture* 
but  does  not  include  software  consulting*  Cost  reduction 
through  Quantity  is  available* 


Fred  Wilke 

TI-59  Custom  Module  Division 

Texas  Instruments 

Lubbock r    TX 

(806)  741-3240 


Robert  Kruger 

Horizons  Technology?  Inc. 

7S30  Clairemont  Mesa  Boulevard 

San  Diegor  CA  92111 

(714)  292-8331 
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